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Motivated by experiments on Josephson junction arrays, and cold atoms in an optical lattice in a 
synthetic magnetic field, we study the "fully frustrated" Bose-Hubbard (FFBH) model with half a 
magnetic flux quantum per plaquette. We obtain the phase diagram of this model on a 2-leg ladder 
at integer filling via the density matrix renormalization group approach, complemented by Monte 
Carlo simulations on an effective classical XY model. The ground state at intermediate correlations 
is consistently shown to be a chiral Mott insulator (CMI) with a gap to all excitations and staggered 
loop currents which spontaneously break time reversal symmetry. We characterize the CMI state as 
a vortex supersolid or an indirect exciton condensate, and discuss various experimental implications. 



The simplest model to understand strongly correlated 
bosons is the Bose-Hubbard (BH) model [l[ which de- 
scribes bosons hopping on a lattice and interacting via a 
local repulsive interaction. With increasing repulsion, at 
integer filling, its ground state undergoes a superfluid to 
Mott insulator quantum phase transition which has been 
studied using ultracold atoms in an optical lattice Q. 

Remarkably, recent experiments have used two-photon 
Raman transitions to create a uniform or staggered "syn- 
thetic magnetic field" for neutral atoms Q, permitting 
one to access large magnetic fields for lattice bosons. The 
multiple degenerate minima in the resulting Hofstadtcr 
spectrum can be populated by non-interacting bosons 
in many ways. Repulsive interactions quench this "ki- 
netic frustration" , leading to unconventional superfluids 
[3-0], or quantum Hall liquids @- Tuning the sign of 
the atom hopping amplitude or populating higher bands 
also leads to such frustrated bosonic fluids [J] . These de- 
velopments motivate us to study the interplay of strong 
correlations and frustration in the fully frustrated Bose- 
Hubbard (FFBH), with half a "magnetic flux" quantum 
per plaquette [a-tZl . At large integer filling, the FFBH is 
also the simplest quantum variant of the classical fully 
frustrated XY (FFXY) model % G3 of Josephson junc- 
tion arrays (JJAs) (Tlj . 

Here, we obtain the phase diagram shown in Fig. Q] of 
the FFBH model at integer filling on a 2-leg ladder us- 
ing the density matrix renormalization group (DMRG) 
method [l2[ and Monte Carlo (MC) simulations. Our 
key result is that the ground state of the FFBH and 
quantum FFXY models at intermediate Hubbard repul- 
sion is a chiral Mott Insulator (CMI). The CMI has a 
nonzero charge gap, and simultaneously supports stag- 
gered loop currents that spontaneously break time rever- 
sal symmetry. With increasing repulsion, the CMI under- 
goes an Ising transition into an ordinary Mott insulator 
(MI) where the loop currents vanish. Weakening the re- 
pulsion leads to a Berezinskii-Kosterlitz-Thoulcss (BKT) 
fl3| transition out of the CMI into a previously studied 
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FIG. 1: (Color online) (A) Phase diagram of the effective clas- 
sical model Hxy, with J T — Ju, obtained via MC simulations 
(see text for details). (B) Phase diagram of the FFBH model 
in Eqn. [T] obtained using DMRG. Both models exhibit a chi- 
ral Mott insulator (CMI) state sandwiched between a chiral 
superfluid (CSF) and an ordinary Mott insulator (MI). (1/J T 
in the XY model ~ sfujl in the FFBH model. [H) 

chiral superfluid (CSF) phase [l4| which retains current 
order. We show that the CMI may be viewed as a vortex 
supersolid or an exciton condensate, and discuss the loop 
current, the charge gap, and the momentum distribution 
across the phase diagram. 

Fully Frustrated Bose-Hubbard Ladder. — The Hamil- 
tonian of the FFBH model on a 2-leg ladder is 

H = -t^{ala x+1 + ai +1 a x )+t^{bib x+l +bi +1 b x ) 

X X 
X X 

where a and b label the two legs of the ladder (see Fig. [2j) , 
t± couples the two legs, and U is the local boson repul- 
sion. The opposite signs of the hopping amplitude (it) 
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FIG. 2: (Color online) (A) Dispersion of the FFBH model at 
U — 0, with two degenerate minima in the low energy a-band. 
Interactions force an equal number of bosons (on average) 
to condense into each of the two minima. (B) Alternating 
pattern of plaquette currents in the presence of chiral order. 



on the two legs leads to an Aharonov-Bohm phase of it 



for a boson hopping around an elementary plaquette [21 



For [7 = 0, the boson dispersion (in Fig. [5] (A)) exhibits 
two bands with the lowest (a) band having degenerate 
minima at momenta fc = 0, ir. This leads to a large degen- 
eracy of many-body ground states — the ground state for 
N bosons corresponds to having N\ bosons in one mini- 
mum and (N—Ni) in the other for any Ni < N — which 
is broken by the repulsion. . The minimum at k = ( k = it) 
has a wavefunction that mainly resides on leg-a (leg-&). 
Since the Hubbard repulsion favors a uniform density, it 
prefers an equal number of bosons at k = 0, ir. A mean 
field Bose condensed state thus takes the form 
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Here if is the U(l) condensate phase, 9 is the relative 
phase between the two modes, and a\ n creates quasi- 
particles at fc = 0, ir. 

For small U, Hartree theory 0, [22j shows 9 = ±ir/2, 
while p has (nonuniversal) power law order. This Lut- 
tinger liquid is the CSF - it supports the long-range stag- 
gered current pattern in Fig. 1(B). The two signs of 9 
correspond to patterns related by time-reversal or unit 
lattice translation. For very large U , both 9 and tp are 
disordered, leading to an ordinary MI which respects all 
the symmetries of H. Remarkably, for intermediate U , 
we find that <p is disordered leading to loss of superflu- 
idity, while 9 is pinned at ±7r/2, spontaneously breaking 
(Ising) time reversal symmetry. This fully gapped inter- 
mediate state is the CMI. This goes beyond mean field 
theory which predicts a direct CSF-MI transition [22j ■ 

Physical pictures for the CMI. — The CSF, with stag- 
gered currents depicted in Fig. [5] (B), is best viewed as 
a vortex crystal where vortices and antivortices are nu- 
cleated by the presence of frustration, and locked into an 
'antifcrromagnctic' pattern due to the intervortex repul- 
sion. At large U, this crystal melts and the vortices com- 
pletely delocalize - this vortex superfluid is well known 
to be simply a dual description of the ordinary MI [27j |. 
However if a small number of defect vortices in the vortex 
crystal delocalize and condense, they kill superfluidity 



but preserve the background vortex crystallinity. This 
vortex supersolid is the dual description of the CMI. 

A different but equivalent picture emerges if we start 
from the usual MI at large U which supports charge 
gapped particle and hole excitations (adding or remov- 
ing bosons). These excitations have degenerate disper- 
sion minima at k = 0,7r as in Fig. 1(A), similar to the 
original noninteracting bosons. Decreasing U decreases 
the MI charge gap. If the charge gap vanishes, the re- 
sulting gapless particles and holes at k = 0, ir could yield 
a Bose condensed (or power-law) superfluid. However, 
a precursor phase emerges from first condensing a neu- 
tral indirect exciton, composed of a particle and a hole 
at different momenta (fc = and k = ir), while the parti- 
cles and holes are still gapped. The CMI is precisely this 
intervening 'exciton condensate' (22|. 

Effective bilayer XY model. — To quantitatively flesh 
out the phase diagram described above, we first study the 
FFBH model at large fillings, where it is equivalent to a 
quantum FFXY model used to describe JJAs of charge 
2e Cooper pairs with an Aharonov-Bohm flux of hc/Ae 
per plaquette. The quantum FFXY Hamiltonian in turn 
maps on to an effective classical model on a 'space-time 
lattice' leading to a classical 2D bilayer square lattice 
model [HI H X y = - J2u J s cos (ip t - ip l+ s) , where ipi 
are the boson phases, and (i, i+S) denote nearest neigh- 
bour sites along 5. The couplings Js take on values ±J\\ 
on the two legs, J± on the rungs linking the two lay- 
ers, and J T in the imaginary time direction [22}. (We 
choose the 'time step' in the imaginary time direction to 
set Jii = J T [22j].) Phase ordering leads to a superfluid, 
while the fully paramagnetic phase of Hxy is the ordi- 
nary MI. Based on small system studies of Hxy [HI, it 
has been argued that the isotropic case Jj_ = Jy exhibits 
a single phase transition with novel exponents, while the 
hig hly anisotropic case harbors two separate transitions 
[231 ]. Here we use extensive MC simulations, on LxLx2 
bilayers with L = 16-64, to obtain the phase diagram 
shown in Fig. 1(A). We find three phases: the CSF, the 
regular MI, and an intervening CMI for a wide range of 
Jj_ including the isotropic point J±_ = Jp . We show that 
CSF-CMI and CMI-MI phase transitions are BKT and 
2D Ising transitions respectively. 

Fig. 3 shows the MC data for J± = l. Similar data was 
also obtained for various J±/J T . Fig. 3(A) shows that the 
helicity modulus T (related to the superfluid density) has 
an increasingly abrupt change with 1/ J T for increasing L, 
indicative of a jump as at a BKT transition. If the transi- 
tion out of the CSF is indeed a BKT transition, V can be 

fit to the finite size scaling form T(L) = A (^1+ 2 (i ^l+c) ) 
(with fit parameters A,C) right at the transition point, 
with A taking on the universal value of 2/7T, while C is a 
non-universal constant 

QUEl. Fitting r(L) to this form, 
we find that the error to this fit shows a sharp minimum 
[13, El at a certain 1/J T (FigfS] inset), with A w 2/tt 
at this dip. This not only allows us to precisely locate 
the transition out of the CSF state, but also confirms its 





FIG. 3: (Color online) (A) Helicity modulus T versus 1/J T 
for different system sizes for Jj_ = 1. (A-Inset) RMS error 
of fit to the BKT finite size scaling form of F shows a deep 
minimum [S^] at the transition, at 1/J T = 0.887(1), and yields 
a jump Ar« 0.637, close to the BKT value 2/tt. (B) Binder 
cumulants for the staggered current versus 1/ J T (for different 
L for J± = 1) intersecting at a continuous transition at 1/ J r = 
0.981(4). (B-inset) Critical susceptibility versus L gives the 
ratio of critical exponents ^/v~ 1.72, very close to 2D Ising 
value 7/1/ = 7/4. Error bars are smaller than the symbol sizes. 



FIG. 4: (Color online) (A) DMRG results for n(k = 0)L" 3/4 
versus U/t, for the FFBH Hamiltonian in Eqn. [T]with t± = 
t and various L. The crossing of these curves at U c i/t ~ 
3.98(1) yields the CMI-MI transition (see text). Inset shows 
the onset of the charge gap at U c i . (B) Rung current structure 
factor Sj(ir)L 20/ " versus U/t at t± = 1. The intersection 
point yields the CMI-MI Ising transition at U C 2 « 4.08(l)t. 
Inset shows Sj{i:)L 2 ^ v versus 5L 1/v with 5 = (U - U c2 )/t, 
for different U/t, leading to a scaling collapse for 2D Ising 
exponents v = 1 and ft = 1/8. 



BKT nature. 

To check for staggered loop currents, we compute the 
Binder cumulant Bl = (l — (TO 4 )i/3(m 2 )|), for the or- 
der parameter m = J2i T i~^T Jiri where Ji tT is the 
current around a spatial plaquette. For small 1/J r , we 
find Bl — > 2/3 indicating long range current order, while 
Bl — > for large 1/J T indicating absence of loop cur- 
rents. Fig. EJB) shows the transition point where the 
current order vanishes as seen from the crossing of Bl 
curves [24| for different L. Remarkably, we find that loop 
current order persists into the regime where the super- 
fluid order is absent, revealing an intermediate insulating 
phase with staggered loop currents - this is the CMI. 

For J±/J = 1, we find the BKT transition occurring 
at 1/J T = 0.887(1) while the current order vanishes at 
the Ising transition which is located at 1/J T = 0.981(4), 
where the error bars on the transition point are estimated 
from the error in the location of the dip in the inset of 
FigE{A) and the error in the crossing point in Fig(3jB), 
both of which yield the limiting thermodynamic values 
for the transition points. This establishes that the phase 
diagram supports three phases: CSF, CMI, and MI. A 
similar analysis for different values of J± allows us to 
obtain the phase diagram in Fig. [HA). 

We have already seen that the transition out of the 
CSF, i.e., the CSF-CMI transition, is of the BKT type. 
The scaling of the divergent susceptibility peak Xait(L) 
for current order (Fig. I3JB) inset) shows that the CMI- 
MI critical point is a 2D Ising transition. Such consec- 
utive, closely spaced, BKT-Ising thermal transitions are 
also observed in the classical 2D FFXY model [X0| . al- 
though its Hamiltonian is quite distinct from Hxy, and 
the chiral order in the classical model corresponds to hav- 
ing in-plane currents rather than interlayer currents as in 
our bilayer model. Such consecutive transitions are also 



found in spinor condensates [25[. 

DMRG study. - We next study the FFBH ladder 
model in Eq. {TJ) at a filling of one boson per site us- 
ing the finite size DMRG (FS-DMRG) method 0. (We 
set t = 1 here.) As noted previously @, [H[, the boson 
momentum distribution n(k) in the presence of 7r-flux 
exhibits two peaks; for our gauge choice, these peaks 
arc located at k = 0, tt. In the CSF state, which is a 
Luttinger liquid [26j on the ladder, we have a singular 
momentum distribution n(k — > 0) ~ \k\~( 1 ~ K / 2 \ with 
K > being an interaction dependent Luttinger param- 
eter [13. Similarly, n(k -> tt) ~ \k - -k\-^- k / 2 \ Let 
U c i denote the location of the transition out of the CSF 
into an insulator. If this transition is of the BKT type, 
as shown from our XY model study, the exponent K 
should take on a universal value K c = l/2 at U c \. A 
plot of n(k = 0)L~ 3 / 4 for different L should thus show a 
crossing point at the transition out of the CSF, as seen 
at U cl w 3.98(1) in Fig. g^A) for t ± = 1. Remarkably, 
Fig. 4(A) (inset) shows that the charge gap also becomes 
nonzero for U > U c \, coinciding with the point where 
A' =1/2, confirming that the CSF-to-insulator transition 
is a BKT transition. This leads to the phase boundary 
of the CSF state shown in Fig. 2(B). 

The staggered current order parameter can be ob- 
tained from the rung-current structure factor Sj(k) — 
T^Y, x ,x' elk(x ^ X ' ) (jxjx'), with j x = i {a\b x - b\a x ) . 
Sj(k = n) ~ L indicates long range staggered current 
order. Our XY model study informs us that the current 
order disappears at a MI-CMI transition which is in the 
Ising universality class. Wc thus expect Sj (tt) to obey the 
critical scaling form Sj(-K)L^/ V = / (([/ - U c2 ) L x l v ) , 
where U C 2 is the CMI-MI critical point, /(.) is a univer- 
sal scaling function, and j3 = 1/8 and v = 1 are the Ising 
critical exponents. As a result, curves of Sj(Tr)L 2 ^» for 
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different L arc expected to intersect at the MI-CMI crit- 
ical point U c2 . This crossing, as seen at U c2 ~ 4.08(1) 
for t±_ = 1 from Fig. HI allows us to carefully locate the 
CMI-MI phase transition. As seen in Fig 2] (inset), plot- 
ting Sjin)^/" as a function of (U - U c <x)lM v shows a 
complete data collapse for U C 2 = 4.08. Similar to our dis- 
cussion for the computations on the XY model, our anal- 
ysis of these crossing points in the FFBH model yields the 
limiting thermodynamic values of the transition points, 
and the error bars are estimated from examining the er- 
rors in these crossing points. Such an analysis, carried 
out for a range of values of t±/t, allows us to map out the 
MI-CMI phase boundary in Fig.QjB); we find U c2 > U cl , 
again consistent with an intermediate CMI state. 

Discussion. — Our computations on the FFBH model 
at unit filling and the XY model (which describes the 
FFBH model at large integer filling), suggest that the 
CMI appears near the tip of the Mott lobes at all boson 
fillin gs o n the ladder. We have generalized the work of 
Ref. [28( to obtain a long-range Jastrow correlated wave- 
function which captures all the essential correlations of 
this CMI state on the ladder [22| ■ Since the CMI is com- 
pletely gapped, with not just a charge gap but also an 
"Ising" gap to charge-neutral excitations, it will be sta- 
ble in a 2D system of weakly coupled FFBH ladders. 

The CSF and CMI states are bosonic an alog s of stag- 
gered current metallic [13] and insulating jl9j | states of 



fcrmions in models of cuprate superconductors. The CSF 
and CMI also find analogs in insulating magnets: param- 
agnetic gapless [HI or spin-gapped [Icj phases with long 
range vector chiral order. 

The CMI may be realized in a Josephson junction lad- 
der at a magnetic field of hc/Ae flux per plaquette fill ) . 
where it would appear as an insulator in transport mea- 
surements. With a Josephson coupling ~ IK, we esti- 
mate that the spontaneous loop currents could produce 
staggered magnetic fields ~ InT for arrays with lattice 
parameter lOum, which could be measured using SQUID 
microscopy [29( . Ultracold bosonic atoms in the presence 
of a (uniform or staggered) synthetic 7r-flux Q are candi- 
dates to realize the CMI. The signature of the flux would 
appear as twin peaks in the atom momentum distribu- 
tion: the peaks would be sharp in the CSF but broad in 
the CMI and MI. Re-interfering the k = and k = ir 
peaks obtained in time of flight via Bragg pulses [3(| 
could test for the persistence of intermode coherence (the 
phase 9 = ±tt/2) in the CMI, and distinguish it from the 
MI. Jaynes-Cummings lattices in a "magnetic field" [3l| . 
could also be used to simulate a polariton FFBH model. 
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